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Abstract 

] The evolution of magnetic fields frozen to a perfectly conducting plasma fluid around a Kerr 

O 

black hole is investigated. We focus on the plunging region between the black hole horizon and 

(N ' 

q ■ the marginally stable circular orbit in the equatorial plane, where the centrifugal force is unable to 

D . 

stably balance the gravitational force. Adopting the kinematic approximation where the dynamical 
effects of magnetic fields on the fluid motion are ignored, we exactly solve Maxwell's equations 
with the assumptions that the geodesic motion of the fluid is stationary and axisymmetric, the 
magnetic field has only radial and azimuthal components and depends only on time and radial 
coordinates. We show that the stationary state of the magnetic field in the plunging region is 
uniquely determined by the boundary conditions at the marginally stable circular orbit. If the 



> 

in 

(N 



magnetic field at the marginally stable circular orbit is in a stationary state, the magnetic field 



Or 



in the plunging region will quickly settle into a stationary state if it is not so initially, in a time 



determined by the dynamical time scale in the plunging region. The radial component of the 
magnetic field at the marginally stable circular orbit is more important than the toroidal component 
in determining the structure and evolution of the magnetic field in the plunging region. Even if 

•rH . 

at the marginally stable circular orbit the toroidal component is zero, in the plunging region a 
toroidal component is quickly generated from the radial component by the shear motion of the 
fluid. Finally, we discuss the dynamical effects of magnetic fields on the motion of the fluid in the 
plunging region. We show that the dynamical effects of magnetic fields are unimportant in the 
plunging region if they are negligible on the marginally stable circular orbit. This supports the 
"no-torque inner boundary condition" of thin disks, contrary to the claim in the recent literature. 
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I. INTRODUCTION 



It is widely believed that accretion disks around Kerr black holes exist in many astro- 
physical environments, ranging from active galactic nuclei to some stellar binary systems 
. People usually assume that the inner boundary of a thin Keplerian disk around a Kerr 
black hole is located at the marginally stable circular orbit, inside which the centrifugal force 
is unable to stably balance the gravity of the central black hole |3], |J . In the disk region, 
particles presumably move on nearly circular orbits with a very small inward radial velocity 
superposed on the circular motion, the gravity of the central black hole is approximately 
balanced by the centrifugal force. As disk particles reach the marginally stable circular or- 
bit, the gravity of the central black hole becomes to dominate over the centrifugal force and 
the particles begin to nearly free-fall inwardly. The motion of fluid particles in the plunging 
region quickly becomes supersonic then the particles loose causal contact with the disk, as a 
result the torque at the inner boundary of the disk is approximately zero |5], ||, |7|, || |9|, and 
references therein]. This is usually called the "no-torque inner boundary condition" of thin 
accretion disks. 

Some recent studies on the magnetohydrodynamics (MHD) of accretion disks have chal- 
lenged the "no-torque inner boundary condition" . Magnetic fields have been demonstrated to 
be the most favorable agent for the viscous torque in an accretion disk transporting angular 
momentum outward [ID, and references therein]. By considering the evolution of magnetic 
fields in the plunging region, Krolik pointed out that magnetic fields can become dynam- 
ically important in the plunging region even though they are not so on the marginally stable 
circular orbit, and argued that the plunging material might exert a torque to the disk at 
the marginally stable circular orbit. With a simplified model, Gammie fl2| solved Maxwell's 
equations in the plunging region and estimated the torque on the marginally stable circular 
orbit. He demonstrated that the torque can be quite large and thus the radiation efficiency 
of the disk can be significantly larger than that for a standard accretion disk where the 
torque at the inner boundary is zero. Furthermore, Agol and Krolik |13| have investigated 
how a non-zero torque at the inner boundary affects the radiation efficiency of a disk. 

Numerical simulations of MHD disks || [14], [15|, 16|, [17| have greatly improved our under- 
standing of disk accretion processes. These simulations show that the magneto-rotational 
instability effectively operates inside the disk and leads to accretion, though the accretion 
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picture is much more complicated than that assumed in the standard theory of accretion 
disks. Generally, the disk accretion is non-axisymmetric and strongly time-dependent. It 
is also found that, as disk material gets into the plunging region, the magnetic stress at 
the marginally stable circular orbit does not vanish but smoothly extends into the plunging 
region [TjJ [IB], [T7|], though the effect is significantly reduced as the thickness of the disk 
goes down || IS|. Furthermore, the specific angular momentum of particles in the plunging 
region does not remain constant, which implies that the magnetic field may be dynami- 



cally important in the plunging region [15, 16, 17]. All these results are fascinating and 
encouraging. Unfortunately, due to the limitation in space resolution and time integration, 
stationary and geometrically thin accretion disks are not accessible to the current 2-D and 
3-D simulations. So it remains unclear how much insights we can get for stationary and 
geometrically thin accretion disks from these simulations ||. 

Instead of small-scale and tangled magnetic fields in an accretion disk transporting an- 
gular momentum within the disk, a large-scale and ordered magnetic field connecting a 
black hole to its disk may also exist and play important roles in transportation of angular 



momentum and energy between the black hole and the disk [18, 19, 20, 21, B3]. Recent 



XMM-Newton observations of some Seyfert galaxies and Galactic black hole candidates pro- 
vide possible evidences for such a magnetic connection between a black hole and its disk 



All these promote the importance of studying the evolution and the dynamical effects of 
magnetic fields around a Kerr black hole. 

In this paper, we use a simple model to study the evolution of magnetic fields in the 
plunging region around a Kerr black hole. We assume that around the black hole the 
spacetime metric is given by the Kerr metric; in the plunging region, which starts at the 
marginally stable circular orbit and ends at the horizon of the black hole, a stationary and 
axisymmetric plasma fluid flows inward along timelike geodesies in a small neighborhood of 
the equatorial plane. The plasma is perfectly conducting and a weak magnetic field is frozen 
to the plasma. The magnetic field and the velocity field have two components: radial and 
azimuthal. We will solve the two-dimensional Maxwell's equations where the magnetic field 
depends on two variables: time and radius, and investigate the evolution of the magnetic 
field. This model is similar to that studied by Gammie \T2\, but here we include the time 
variable. Furthermore, we ignore the back-reaction of the magnetic field on the motion of 
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the plasma fluid to make the model self-consistent, since if the dynamical effects of the 
magnetic field are important the strong electromagnetic force will make the fluid expand in 
the vertical direction. The ignorance of the back-reaction of the magnetic field will allow 
us to to analytically study the evolution of the magnetic field, but it will prevent us from 
quantitatively studying the dynamical effects of the magnetic field. However, we believe 
that the essential features of the evolution of the magnetic field in the plunging region is not 
sensitive to the details of the dynamical effects. To check the self-consistency of the model, 
we will estimate the dynamical effects of the magnetic field by considering the back-reaction 
of the magnetic field on the fluid motion. We will also discuss the self-consistent solutions to 
the coupled Maxwell and dynamical equations, and look for implications to the "no-torque 
inner boundary condition" . 

The paper is organized as follows: In Sec. |I| we write down the Maxwell's equations for 
an ideal MHD plasma. In Sec. [1 1 J we write down the general forms of the magnetic field 



and the velocity field of the plasma fluid around a Kerr black hole. In Sec. [TV] we solve 
Maxwell's equations with the approximations outlined above. In Sec. |V], with the solutions 
obtained in Sec. we study the evolution of the magnetic field in the plunging region. 
In Sec. [VT|, we check if magnetic fields can become dynamically important in the plunging 
region. In Sec. |V1J| we draw our conclusions. In the Appendix we solve Maxwell's equations 
in the disk region, since we want our solutions in the plunging region to be continuously 
joined to the solutions in the disk region. 

Throughout the paper we use the geometrized units G = c = 1 and the Boyer-Lindquist 



coordinates (i, r, 8, <p) p6| £7] ■> except in the Appendix where cylindrical coordinates are 
used. 

II. MAXWELL'S EQUATIONS FOR AN IDEAL MHD FLUID IN A CURVED 
SPACETIME 

In a curved spacetime, Maxwell's equations are 

V [a F bc] = , (lb) 
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where F a b = —Fb a is the electromagnetic field tensor, J a is the current density 4-vector of 
electric charge. In Eqs. (|la|) and (TTB]), V a denotes the covariant derivative operator that is 



compatible with metric: V ' a 9bc = 0; the square brackets "[ ]" denote antisymmetrization of 
a tensor. 

For an ideal MHD plasma fluid whose electric resistivity is zero, the electric field in the 
comoving frame is zero, i.e. 

E a = F ab u b = , (2) 

where u a is the 4-velocity of the fluid. In other words, the magnetic field is frozen to the 
plasma fluid. Then, the electromagnetic field tensor F a b can be written as 

F a b = —£ a bcdB c u d , (3) 

where B a is the magnetic field measured by an observer comoving with the fluid (i.e., having 
a 4-velocity u a ), and e a b c d is the totally antisymmetric tensor of the volume element that is 
associated with the metric g a b- By definition, the magnetic field B a satisfies 

B a u a = . (4) 

The corresponding electromagnetic stress-energy tensor is 

T ™ = ~ \9 ab FcdF cd ^j = ^B\ a u b + ^B 2 g ab - ±-B a B b , (5) 

where B 2 = B a B a . 

In the case of MHD, the electric current density J a is unknown, but it is only defined 



by Eq. (la). The fundamental variables are the electric field E a and the magnetic field B a . 
Maxwell's equations are then reduced to Eq. (|TbD . Since in the comoving frame of an ideal 
MHD flow E a = and B a is related to the electromagnetic tensor F a b by Eq. (R|), the dual 



of Eq. (TLB gives p| 



V a (u^) = . (6) 



For an ideal MHD fluid Eqs. ( pJbf) and @ are equivalent. 
Eq. (H) can be expanded as 



u a V a B b = B a V a u b - B b V a u a + u b V a B a . (7) 



The contraction of Eq. ([7]) with u b leads to 

V a B a - B a a a = , (8) 

where a a = u b V b u a is the acceleration of the fluid. In deriving Eq. (||) we have used Eq. (f|) 
and the identity u a u a = — 1. 

Substituting Eq. (g) into Eq. (0), we get 

u a V a B h = B a V a u b - B b V a u a + u b B a a a . (9) 

The tensor V a u b is decomposed as 



V a M b = ^Qh ab + a ab - uj ab - u a a b , (10) 

where h ab = g a b + u a Ub is the space-projection tensor, 9 = h ab V a u b = V a u a is the expansion, 
cr ab = V( a ii&) — Qh ab /3 is the shear tensor, and u; a fe = — V[ a u b ] is the vorticity tensor of the 
fluid. Here the braces "( )" denote symmetrization of a tensor. It is easy to check that 

a ab u b = u ab u b = , a ab h ab = a ab g ab = . (11) 

Substituting Eq. (|10"D into Eq. (|), we obtain 

u a \7 a B b = B a (a ab - u ab ) - ^OB b + B a a a u b , (12) 

which shows that the evolution of the magnetic field is governed by the expansion, the shear, 
the vorticity, and the acceleration of the fluid. The contraction of Eq. (O) with B b gives 



^u a V a B 2 = B a B b a ab -^QB 2 , (13) 

where we have used B a B b u) ab = and Eq. (^). 

We note that with the space-projection tensor h ab , Eq. can be written as 

h ab V a B b = , (14) 

which says that the spatial divergence of the magnetic field is zero. 

In terms of differential forms, the tensor F ab is a closed 2-form, then the Maxwell equation 
( |lb"l) can be written as dF = [|7|, |30|. Using Stokes' theorem and Eq. (H), it can be shown 



that for a perfect conducting fluid the magnetic flux threading any 2-dimensional spatial 
surface S is unchanged as the surface moves with the fluid 

$ B = / F = constant . (15) 
Js 

This is the mathematical formulation for the statement that magnetic field lines are frozen 
to a perfectly conducting fluid. 



III. PARTICLE KINEMATICS AND MAGNETIC FIELDS AROUND A KERR 
BLACK HOLE 



Now let us assume that the background spacetime is the outside of a Kerr black hole of 
mass M and angular momentum Ma, where — M < a < M. In Boyer-Lindquist coordinates, 
the Kerr metric is EI El] 



, 2Mr\ , 9 AMar „ „ , , £ , , Asm 2 9 , 9 

rfs 2 = - ( 1 - — — dt 2 — sin 2 9 dtdcj) + — dr 2 + Zd9 2 + — - — d(f> 2 , (16) 

2j / Z.J Za 2-j 



where 

A = r 2 - 2Mr + a 2 , £ = r 2 + a 2 cos 2 9 , A = (r 2 + a 2 ) 2 - Aa 2 sin 2 5 . (17) 
A Kerr black hole usually has two event horizons: an inner event horizon and an outer event 



horizon, whose radii are given by the two roots of A = |37| • What is relevant in this paper 
is the outer event horizon, so whenever we talk about the "event horizon" we always mean 
the outer event horizon, whose radius is 

r H = M + (M 2 - a 2 ) 1/2 . (18) 

Let us define an orthonormal tetrad attached to an observer comoving with the frame 
dragging of the Kerr black hole, {eg, e*, e^, eg}, by 

„ 1 



d\ a ( d 
dt ■ ; 



V2 (°X e a = — (—Y e«=(*) — (± ) (19) 
sJ [dr ' 2 ~EV 2 [d6 ' 3 ~[A sm9[d^> ' ' ' 



where 



which are respectively the lapse function and the frame dragging angular velocity. As r — » rn, 
we have x an d u ~^ ^h, where 

Q H = (21) 

H 2 Mr n V ; 

is the angular velocity of the event horizon. Then, the 4- velocity of the fluid, u a , can be 

decomposed as 

u a = r(e% + v f el + v s eZ + v$ei) , (22) 



where (v?,vg,v£) are the components of the 3-velocity of the fluid relative to the observer 
comoving with the frame dragging, and F = (1 — v% — f| — ^|) -1 ^ 2 is the Lorentz factor. 
Inserting Eq. ( |T9"D into Eq. (^), we have 



u 



d_ 

X \dt 



'A 
Vf I — 



1/2 



d_ 

dr 



d 



E 1 / 2 V 09 



X \d(j) 



where 



Q = u + xv^ 
is the angular velocity of the fluid. 



.4 



1/2 



sin 9 



(23) 



(24) 



The specific angular momentum of a fluid particle is 



U r . 



TA sin 2 9 



(Sl-w) = Tvi - sin#. 



.4 



1/2 



(25) 



Obviously, L = when Q = uj. Thus, an observer comoving with the frame dragging of a 



Kerr black hole has zero angular momentum plj . The specific energy of a fluid particle is 

d 



E 



Of 



Un 



Tx + Loo. (26) 

In Eq. (p6|), the term Luo represents the coupling between the orbital angular momentum of 
the particle and the frame dragging of the Kerr black hole. When the particle moves on a 



geodesic, L and E defined above are conserved [^g, \Z]\, |29|| . 

The magnetic field B a , which satisfies Eq. (|j), can be decomposed as 

B a = (B f Vf + B §V§ + B^)e a + B f e\ + B § e a 2 + B $ e a 3 , 

from which we have 



B 1 



B a B a = Bj + B\ + B\- (B f .vr + B § v § + 5 



XV 1 



(27) 



(28) 



Inserting Eq. (pT9|) into Eq. (^7]), we have 



B a = —(BfV? + B§v§ + B 



X 



M> 



d 

3t ] 



Ba 



d 



E 1 / 2 V 39 
and, correspondingly 



7777 + 



<P_ 

sin 9 \A 



d_ 

ft 

1/2 / d 



Bf I ^ 



1/2 



d_ 

Or 



(29) 



B n 



1/2 

-x(B r Vf + 5^ + B^vj) dt a + B r ( — ) dr a 



/4\ 1/2 

+^E 1 / 2 ^ a + ^(-j sin# 



(d0 a - u;dt a ) . 



(30) 
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IV. SOLUTIONS OF MAXWELL'S EQUATIONS 



The Maxwell equations that we want to solve are given by Eq. . In terms of coordinate 
components, Eq. (|]) takes a very simple form 

1 %[V=g(^-^B a )]=0 1 (31) 



y/—g dx c 

where x a = (t,r,6,<p), g is the determinant of the metric tensor g Q p, u a = u a dx^ and 
B a = B a dx2 are respectively 

= - r = r^(|) 1/2 , «' = |^, = (32) 

and 

/a\ 1/2 R. 5~ /y\ 1/2 

• B '^' + ^*' (33) 

Because of the constraint Eq. (^), among the four equations of (EH]) only three are indepen- 
dent. 

Since the background spacetime is stationary and axisymmetric, we can look for station- 
ary and axisymmetric solutions with d/dt = 8/ dip = 0. However, since we are interested 
in the time evolution of magnetic fields, we will keep the d/dt terms on magnetic fields but 
adopt that d/d<p = 0. To simplify the problem, we further assume that in a small neighbor- 
hood of the equatorial plane (i.e., \ir/2 — 9\ <C 1), u 6 = B 9 = (i.e., v § = B§ = 0). This 
assumption, which has also been used by Gammie fT2] , ensures that du e /d6 = dB e /d8 = 
on the equatorial plane. We emphasize that, when d/dt ^ 0, this assumption holds only 
if the fluid moves geodesically, which requires that the magnetic fields are weak and their 
dynamical effects are negligible. Otherwise the electromagnetic force will make u 9 and B e 
non-zero except exactly on the equatorial plane. Thus, hereafter we assume that fluid par- 
ticles move on timelike geodesies in the plunging region. This assumption will be justified 
latter. 

Now, let us focus on solutions on the equatorial plane (6 = tt/2). Considering the fact 
that for the Kerr metric \f—g = S sin 9 = r 2 on the equatorial plane, Eq. ( j3~l~|) is reduced to 

y^B? - u^B 1 )] + [r 2 (u r B? - u?B r )] = . (34) 

9 



For (3 — r,t, Eq. (|34|) gives 

-| [r 2 (n r 5* - u'B 7 ')] = 
respectively. The solution of Eq. fl35|) is 



<9r 



(35) 



u l B r - u r B 



t Co 



(36) 



where Cq is a constant. With Eq. (|T|), it can be checked that 27rCoA# is the magnetic flux 
threading the 2-dimensional surface defined by t = constant, r = constant, || — d\ < ^ 
and < 6 < 2n. 



For (3 = 9, Eq. (34) is automatically satisfied since u = in a small neighborhood 

of the equatorial plane everywhere and all the time. 
For $ = (/), Eq. O gives 



at or L 



(37) 



From the constraint Eq. (HI), we have 



£? 9 



(B r u r + E^ut) 



(38) 



Substituting Eq. fl3"B"|) for 5* into Eq. (^), we obtain 



{u l u t + u r u r )B r - 



C u t 



(39) 



Now, substitute Eq. ( p6[ ) for 5* and Eq. (p9|) for .B^ into Eq. (|3T|), we obtain a first order 
partial differential equation 



(40) 



where 



V(t,r) = — (r 2 B r + C u t ). 



(41) 



In deriving Eq. ( H0|) we have used du a /dt = 0. Eq. (Bjj) simply says that \I/ is conserved 
along the worldline of a fluid particle: u a V a ^ = 0. Let us define 



r r = 



XT 



Mi 



(42) 
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which is the coordinate time interval spent by a fluid particle to move from r = ro to r. 
Then Eq. fl4"0|) can be written as 

l + *<*•'> = °- (43) 

The solution of Eq. flUf) is simply 

W = (t - r) , (44) 

i.e. \I/ is a function of t — r. 

Eq. (0) gives a "retarded" solution to Eq. 04"3|) : at any time t the solution at radius r is 
given by the solution at an earlier time t — r(r) at the radius ro- Thus, a variation in the 
magnetic fields at any r will propagate with the fluid motion (Fig. p]). The solution is unique 
if a suitable boundary or initial condition is imposed. For example, if a boundary condition 
is given on t = 0: ty(t = 0, r) = ^i(r), then the solution is ty(t, r) = \?i(t — £). In order for 
the solution to exist for a region specified by t > in the t, r-space, the boundary function 
\I/i(t) must be defined on the whole r-axis: — oo < r < oo. Similarly, if a boundary condition 
is given on r = (i.e., r = ro): ty(t, r = 0) = ^(i), then the solution is \l/(t, r) = ^ 2 (t — r). 
In this case, in order for the solution to exist for a region specified by r > in the t, r-space, 
the boundary function ^^(t) must be defined on the whole t-axis: — oo < t < oo. 

We can also specify the boundary condition in another way 

#(* = 0, r > 0) = ^i(r) , #(t > 0, r = 0) = # 2 (£) . (45) 
Then, the solution of Eq. (|43|) is 

f *i(T-t) , < t < r 
* = < . (46) 

[ * 2 (t-r) , < r < t 

I.e., the value of \& in region < t < r (region I) is determined by the value of \I/ on the 
boundary t = 0, r > 0; the value of \l/ in region < r < t (region II) is determined by 
the value of ^ on the boundary r = 0,t > (see Fig. §). In order for the solutions to be 
smoothly matched on the diagonal line t = r > separating region I and region II, ^i(r) 
and ^2{t) must be smoothly matched at t = r = 0: 



*,«,) = * 2 (0), d * lM 



c/r 
11 



cNf 2 (t) 



r=0 dt 



(47) 

i=0 



Given the solution of \I/ in Eq. (f44|), we can solve B r from Eq. (|41|), then solve .B^ 1 from 
Eq. (p9|). The results are 

B r = -3 (-C wt + ^0) , (48a) 
= J-I-Coutu* + (1 + . (48b) 

Using Eq. fl36|), we obtain 

1 

r 2 ii' 

Note that 5*, £T, and satisfy the constraint Eq. (|38|) , so among the three components 
only two are independent. 

From Eqs. ©, (ggg), (pED , and the fact that u t = u a {d/dt) a = -E and = 



B l = -2-KI + wV)C + • ( 49 ) 



u a {d/d(j)) a = L, we can solve for B? and Si 

C E + . , 

5 < = rA i/ 2 . ( 50a ) 

g (^ + ^)(»-^)r + JL, (50b) 

<p Exrvf Ervf 

where Eqs. ( p5|) and (|26|) have been used to simplify the expression for Si. Since we focus 
on the solutions on the equatorial plane, here and hereafter we set S = r 2 , and use x, w, 
and A to refer their values at 9 = ir/2. Note, in the solutions in Eqs. (|48a ) and ( 48b ) [or, 



equivalently, Eqs. ( |50a| ) and ( p0b| )] all the dependence on time t is contained in the function 

Since we assume that the fluid particles move on geodesies, the specific angular momen- 
tum L and the specific energy E are constants. If \I/ is also a constant, the combination 
CqE + is a constant, which we denote as C\. Then, Eqs. (|50a[ ) and ( |50bQ become 



C 

B * = j^m • < 51a ) 

which are stationary solutions of Maxwell's equations. 

If the black hole is a Schwarzschild black hole (i.e., the specific angular momentum a = 0), 
the stationary solutions are reduced to 

B f = -4 , (52a) 
XT 2 

B * = I^: + I^- (52b) 
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where E = Fx, L = Tx' 1 ^ 2 , and x 2 = I - 2M/r. The non-stationary solutions can be 
obtained by replacing C\ with CqE + ^/L. 



V. EVOLUTION OF MAGNETIC FIELDS IN THE PLUNGING REGION 

We are interested in the evolution of magnetic fields in the plunging region in the equa- 
torial plane between r = r ms , the marginally stable circular orbit, and r = rn, the event 
horizon of the black hole. 

For direct circular orbits (i.e., corotating with L > 0) around a Kerr black hole, the radius 
of the marginally stable orbit is given by [|3T 



M {3 + z 2 - sign(a) [(3 - Zl )(3 + z x + 2z 2 )] 1/2 } 



(53) 



where sign(a) = 1 (0, —1) if a > (= 0, < 0), and 



M 2 ) 

,2 \ 1/2 



ax 1 / 3 / a x 1 / 3 

1 + m) + ( 1 -m) 



^2 = 4+3 



M 2 



(54a) 
(54b) 



The angular velocity of a particle geodesically moving on the marginally stable circular orbit 
is 



J ' ms I r 3 



. M 

l + o -s- 



1/2' 



(55) 



From Eqs. 
circular orbit by 



, we can calculate the circular velocity on the marginally stable 



4>,ms 



A 1/2 
^msXmi 



(^ms ^ms) 



(56) 



the corresponding Lorentz factor by r ms = (1 — ^l^) -1 ^ 2 ) specific angular momentum L ms 
by Eq. (p5[ ), and specific energy S ms by Eq. (^) (setting 6 = -n/2 and r = r ms ). Here and 
hereafter the subscript "ms" represents the values at r = r ms . 

Now let us choose the boundary radius r = r ms , the conserved specific angular momen- 
tum and specific energy to be 



L = L ms (l -6), E = E n 



(57) 
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where < 5 C 1. I.e., we keep the specific energy the same as that on the marginally stable 
circular orbit, but decrease the specific angular momentum by a little amount. Then, we 
can calculate the Lorentz factor at r = r ms , corresponding to the specific energy and the 
specific energy specified by Eq. (|5"7|), by 

r = —{E - Lu ms ) ■ (58) 

Xms 

and the corresponding boundary values of vi and v? at r = r ms by 

Lr ( 1 \ 1/2 

^o = -^' ^o = -(l-^| -p2 j • (59) 

From Eqs. Q25D and (^), we can calculate T, Q, vi, and at any r by 

T = -(E-Luj) , (60a) 

+ (60b) 

V * = A^(E-Lu) ' (6 ° C) 
Y 2 / L 2 r 2 \l 1/2 

= + — )J • < 60d > 

where E and L are given by Eq. fl57|). The parameter r defined by Eq. ( pE2|) can then be 
calculated with 



r(,) = l.^*' (61) 

With the above formulae at hands, we can calculate B?, B^, and B 2 = B 2 + B-, — (BfVf + 
B^v^) 2 at any r by Eqs. ( |50a|) and (|50b|) , giving the constant C and the function ^/(t — r). 
To determine Co and we need to specify the boundary conditions for the magnetic field. 
We will consider stationary solutions and non-stationary solutions separately. 



A. Stationary solutions 

It is straightforward to specify the boundary conditions for stationary solutions (i.e. 
solutions with \l/ = constant). To determine the solutions, we need only specify the values 
of Bf. and B^ at r = r ms : B? and B^ Q . Then, by Eqs. ( |51a| ) and ( |51b| ), we have 



Ci = ^oXmsA^ 2 , (62a) 
^ = B $Q Er ms Vro ~ B m T A]/ 2 {Q - u ms v 2 Q ) , (62b) 
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where G\ = CqE + ^/L, Q is the angular velocity corresponding to 



up O^msX 



^ru' msAms /nn\ 

il = uj mB -\ — 2 . (63) 



inis 



With the G\ and \I/ determined above, we can calculate .Bf, -B^, and B 2 = B 2 + 5| — 
(S^ff + B^v^) 2 at any r by Eqs. ( |51a| ) and (|51b| ). Since Bf and B^ linearly depend on Bfg 
and _B^ , it is sufficient to study the effects of Bf and B^ Q separately. The results for any 
linear combination of B?o and _B^ are simply linear superpositions of the results for Bfo and 
B^ Q separately. 

In Fig. we show the evolution results of B^ with the boundary condition Bi = B^ 
and Bf = at r = r ms , for different spinning status of the black hole and different values 
of 5 that specify the kinetic boundary conditions of the fluid. All quantities are scaled to 
the mass of the black hole so we do not need to specify the value of M. Since Bf does not 
depend on B^ , Bf is always zero. Since C\ = and E = constant, B^ evolves according 
to B^ oc (rvf)~ l . Though in the plunging region r decreases, \v?\ grows faster except at the 
neighbor of r ms . So the evolution of Bi is dominated by the variation of Vf. Thus, as fluid 
particles get into the plunging region, decreases quickly as clearly shown in Fig. ||. This 
radial expansion effect is not sensitive to the spin of the black hole, but very sensitive to the 
value of 5 (or, effectively, the initial value of Vf). As 5 decreases (i.e., \v? \ oc S 1 ^ 2 decreases), 
the fluid expands more as it gets into the plunging region, so the value of B^ decreases more. 

In Fig. |4] we show the evolution results of Bf (dashed lines) and B^ (solid lines) with 
the boundary condition Bf = Bf and B^ = at r = r ms . Though at r = r ms we have 
B^ = 0, in the plunging region B^ becomes nonzero since B^ depends on both Bfo and B^ Q 
[Eqs. ( |51b| ), ( |S2a| ) and ( 62b )]. This is the manifestation that the shear motion of the fluid 



in the plunging region generates B^ from Bf. The radial component of the magnetic field, 
Bf, increases gradually as the fluid enters the plunging region, according to Bf cx (rA 1//2 ) _1 , 
and blows up on the black hole horizon where A = 0. The shear motion of the fluid does 
not amplify Bf, which echos with the fact that Bf is always zero if it is zero at r = r ms 
(Fig. |H). Since Bf does not depend on Vf, there is only one dashed line in each panel of 
Fig. f|. In comparison, the toroidal component, B^, increases more quickly in the transition 
region, since the shear motion of the fluid magnifies B^. This is more prominent for small 
values of 5 (i.e. small \vfo\), since Bi oc Vf 1 [Eq. (|51b|) 1 and Vf is close to Vfo as the particles 



just leave the marginally stable circular orbit. For small values of 5, B^ increases sharply 
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as the fluid just gets into the plunging region, then decreases a little bit due to the radial 
expansion of the fluid. Unlike B? , is always finite on the black hole horizon. We see that, 
Bf. and Bi evolves in very different ways. 

From Fig. |] we also see that the evolution of the magnetic field in the plunging region 
depends on the spin of the black hole, though not very sensitively. Interestingly, the toroidal 
component and the poloidal component depend on the spin of the black hole in an opposite 
way. As the dimensionless spin parameter a/M increases from zero to positive values, the 
shear amplification effect on the toroidal component of the magnetic field increases (except 
for the case of 5 = 1CT 2 for which the amplification effect is not prominent), while the 
amplification of the radial component caused by the compression in the azimuthal direction 
(i.e., the decrease in radius r) decreases. But, if a/M decreases from zero to negative values, 
the shear amplification effect on the toroidal component of the magnetic field decreases, while 
the compression amplification of the radial component increases. The opposite dependence 
for positive and negative spins is probably due to the different coordinate distances from 
the marginally stable circular orbit to the black hole horizon for black holes of positive and 
negative spins. 

In Fig. [| we show the evolution of B 2 [defined by Eq. (p5)j with the same boundary 
conditions in Fig. [|. As the fluid just enters the plunging region (near the right ends of 
curves), B 2 sharply increases due to the small values of \v?\ there, which is the manifestation 
of amplification effect caused by the shear rotation of the fluid. After that, i.e., after the 
fluid obtains a large radial velocity (the dashed lines in the figure), the amplification effect 
is reduced but the expansion effect becomes prominent [see Eq. ([13])]. On the horizon of 
the black hole B 2 is always finite, so the boundary conditions on the horizon is satisfied |32|, 
and references therein]. 



B. Non-stationary solutions 

To specify the boundary conditions for non-stationary solutions is a little bit complicated. 
As discussed earlier, to determine the solutions in a region with t > and r > (i.e. 
?"h < r < r ms ) in the t, r-space, we need to specify the boundary conditions on the axes 
(t = 0,r > 0) and (t > 0,r = 0), i.e. specify * x (r) and tf a (t) [see Eq. (||)]. 

As an example, let us assume that — constant, and Bi = on the axis (t = 0, r > 
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0). I.e., the solution on the boundary r = is stationary for t > 0, and the 0-component 
of the magnetic field in the plunging region is zero at t — 0. Equivalently, we specify the 
boundary conditions as follows: 

On r = r ms , t > : B f = B f0 , B^ = 0; (64a) 

On t = 0, r H < r < r ms : £7 = . (64b) 



Then, if we apply Eqs. (|50a|) and ( |50b|) to the boundary at r = r ms and t > 0, we can solve 
for Co and \I/2- The results are 

C = r [(ffo - ^ ms ^o)^ r o + Xms] , (65a) 
^ 2 = -r 0J B, ^/ s 2 (fi - uj ms v 2 f0 ) . (65b) 



Applying Eq. ( p0b| ) to the boundary at t = and rn < r < r ms , and substituting Eq. ( |65a| ) 
for Co, we can solve for 

*i(t) = *i[r(T)] , * lW .|£^«, («,) 

where r(r) is given by the inverse of r(r). The values of Bf on t = and rn < r < r ms is 
then determined by Eq. ( |50a|) . 



With the and ^2 determined above, we can obtain \1/ by Eq. (46). This, together with 



the C given by Eq. (|65a|) , allows us to calculate Bf and Bi with Eqs. ( |50a|) and ( |50b|) for 



any radius in the plunging region at any time t > 0. 

It is easy to check that ^(r = 0) = ^2, so the solutions are continuously matched on 
the line t = r. Since 

(i*i(r) d\^i(r) dr d^^r) x^fA 1/2 



dr dr dr dr 



(67) 



the solutions are smoothly matched on the line t = r only in the limit Vfo = since then 
we have d^i/dr = = —d^/2/dt at r = t — 0. If Vfo is not zero but small (|ffo| 1), the 
solutions are approximately smoothly matched on the line t = r if |rf^i(r)/dr| is not large 
at v — r ms . 

In Figs. P - |] we show the results for the non- stationary evolution of magnetic fields. 
Each figure corresponds to a different spinning state of the black hole. In Fig. |6] a/M = 0, in 
Fig. [7| a/M = 0.99, while in Fig. |8| a/M = —0.9. The boundary conditions for the magnetic 
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field are given by Eqs. (|64a ) and ( |64b|) . The kinetic boundary conditions of the fluid are 
given by Eq. (0), here we assume 5 = 10~ 3 . In each figure, each panel corresponds to 
the state of the magnetic field at a particular moment. Especially, the first (left and up) 
panel shows the initial (t = 0) state of the magnetic field [i.e. Eqs. fl64ap and ( J64bp ]. In 
each panel, the thick dashed line shows B?, the thick solid line shows Bi, and the thin lines 
show the corresponding stationary solutions. Each curve starts from the marginally stable 
circular orbit (right end) and ends on the horizon of the black hole (left end). 

From these figures we see that, though initially the magnetic field is deviated from the 
stationary state, it evolves to and finally saturates at the stationary state as time goes on. 
Since we hold the magnetic field on the marginally stable circular orbit in a stationary state 
for t > [Eq. (|64a )] and the fluid moves from the marginally stable circular orbit toward 
the horizon of the black hole, the stationary state propagates from the marginally stable 
circular orbit toward the horizon of the black hole (i.e., from right to left in the figures). 
This is clear seen in Figs. ^| - [8]: the magnetic field at a radius closer to the right end of 
the curve gets into the stationary state earlier. In each panel, with a black dot we show the 
position of a particle that is at the marginally stable radius at t = 0, which clearly shows 
the propagation of the state of the magnetic field with the motion of the fluid. This suggests 
that the stationary state of the magnetic field in the plunging region is uniquely determined 
by the boundary conditions at the marginally stable circular orbit. 

The figures also show the dependence of the evolution of the magnetic field on the spin 
of the black hole. For a black hole with a negative a a longer coordinate time is needed for 
most of the fluid in the plunging region to settle into the stationary state. This is apparently 
caused by the fact that a black hole with a negative a has a larger coordinate radial distance 
from the marginally stable circular orbit to the horizon. 



VI. CAN MAGNETIC FIELDS BECOME DYNAMICALLY IMPORTANT IN 
THE PLUNGING REGION? 

In previous sections we have studied the evolution of magnetic fields in the plunging region 
and seen that magnetic fields can be amplified by the convergence (the radial magnetic field) 
and shear motion (the toroidal magnetic field) of the fluid. Thus, we can ask if magnetic fields 
can become dynamically important in the plunging region assuming that their dynamical 
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effects are negligible at the marginally stable circular orbit. We try to answer this question 
in this section. 



A. Dynamical effects of the magnetic field 

Assuming that in the plunging region the gas pressure of the fluid is negligible, then the 
total stress-energy tensor of the fluid is 

T ab = p m u a u b + T£ b A , (68) 

where p m is the mass-energy density of the fluid matter measured by an observer comoving 
with the fluid, is the stress-energy tensor of the electromagnetic field as given by Eq. (|5|). 
Then, the equation of motion of the fluid is given by 

V a T ab = Pm u a V a u b + u b V a (p m u a ) + V a T{? M = . (69) 

When Maxwell's equations are satisfied we have V a T|^ = -F ab J b |2Cj. Then, by Eq. (@), 
we have UbVaT^ = for a magnetic field frozen to the fluid. Therefore, the contraction of 
Ub with Eq. (|69[) leads to the equation of continuity fl38| 

V J£ = , J^ = p m u a , (70) 

where is the flux density vector of mass. The equation of continuity does not explicitly 
contain magnetic field variables. Since (d/dt) a and (d/d<f)) a are Killing vectors, Eq. ( j69|) 
leads to the conservation of angular momentum and energy 

V a J L a = , V a J£ = , (71) 

where 

J ^(^)'^=("- + s B2 ) £ "°-^ B ° (72) 

is the flux density vector of angular momentum, 

% = - (I) V = (p m + ^ Eu« + ^B t B« (73) 
is the flux density vector of energy. 
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Here we look for stationary and axisymmetric solutions in the plunging region. Then, 
with the assumptions adopted in this paper, i.e. d/dt = d/d(f> = 0, u 6 = B 9 = and 
du 6 /d0 = dB e /d6 = 0, in the equatorial plane the equation of continuity is reduced to 

d 



dr 



(r 2 p m u r ) = . 



(74) 



Similarly, the equations of angular momentum and energy conservation are reduced to 



and 



d 
dr 



d 
dr 



—B- I Lu r - —BJV 



An / An 



(75) 



— B 2 Eu r + — B t B 

47T / An 



. 



The solution of Eq. ([[4]) is 



(76) 



(77) 



Anr 2 u r A7rTxA 1 / 2 Vf ' 
where F m = Anr 2 J^ is a constant to be determined by the boundary conditions, which 
measures the mass flux across radius r. Fig. § shows the variation of the mass density 
with radius in the plunging region, assuming that fluid particles move on geodesies. The 
kinetic boundary conditions are given by Eq. (|57|). In each panel (corresponding to a different 
spinning state of the black hole), we show the ratio of p m / p m0 corresponding to four different 
values of 6: 10~ 2 , 10" 3 , 10~ 4 and 10" 5 , where p m0 = p(r = r ms ). We see that, the evolution 
of the mass density of the fluid sensitively depends on the value of 5 (or effectively, the 
value of \vfo\ oc S 1 ^ 2 ). While for not very small 5 (e.g., 5 = 1CT 2 ) the variation of p m is 
not dramatic, for very small 5 (e.g., 5 = 10 -5 ) the variation of p m is dramatic: as the fluid 
gets into the plunging region the mass density drops sharply. This is caused by the sharp 
increase in the ratio Vf/vfQ in the plunging region for small 5. 

From Eqs. (|75j) and ( [76]) we see that the dynamical role played by the magnetic field is 
characterized by the following three dimensionless parameters 



B 2 



Anpv 



V2 



Ba,B 7 



V3 



B t B r 



Anp m Eu r 



(78) 



A r Kp ra Lu r 

The parameter 771 measures the transfer of angular momentum and energy from the fluid 
particles to the magnetic field; 772 and 773 measure the transportation of angular momentum 
and energy from one part to another by the magnetic tension in the rest frame of the fluid. 
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To justify the assumption that the effects of the magnetic field on the dynamics of the fluid 
particles are negligible, we must require that 

7/1 < i , m < 1 , 7/3 « 1 • (79) 

In Eq. (0), p m is given by Eq. Q77p, 5 2 is calculated by Eq. @D (setting ^ = 5^ = 0), 
B^B r and -B t -B r are calculated by 

B^B r = ^B f B $ , B t B r = -*^-B f (xrB?v? + A^Bfil) . (80) 

Near the marginally stable circular orbit, we have \vf/v^\ <C 1, E > LQ, and B t B r « 
—QB ( j > B r . Therefore, from the definitions of 7/2 and 7/3, we have 

7/3 « —7] 2 < 7/2 , (81) 

at r « r ms . 

From Fig. |9| we see that for small 5 the mass density p m drops quickly as the fluid enters 
the plunging region. However, the evolution of B 2 in the plunging region sensitively depends 
on the orientation of the magnetic field on the marginally stable circular orbit. From Fig. |5|, 
the value of B 2 corresponding to an initially radial magnetic field increases as the fluid enters 
the plunging region, because of the convergence of the radial magnetic field and the shear 
amplification of the toroidal magnetic field. If the initial magnetic field is purely toroidal, 
on the other hand, from Fig. |] we see that the magnetic field keeps purely toroidal in the 
plunging region and the dilution effect arising from the drop in the mass density makes 
B^ (and thus B 2 ) decrease in the plunging region. Therefore, in the plunging region the 
variation of rji, 7/ 2 , and 7/3 sensitively depends on the boundary condition of the magnetic 
field on the marginally stable circular orbit. We can imagine that, if on the marginally stable 
circular orbit the magnetic field is purely radial, then in the plunging region B 2 increases so 
the ratio 7/1 also increases as p m decreases. In such a case, 771 may become close to or even 
greater than 1 in the plunging region even if it is 1 on the marginally stable circular orbit, 
then the dynamical effects of the magnetic field become important in the plunging region. 
On the other hand, if on the marginally stable circular orbit the magnetic field is purely 
toroidal, then in the plunging region B 2 decreases which causes 7/1 to decrease also if B 2 
decreases faster than p m . Therefore, to correctly estimate the dynamical effects of magnetic 
fields in the plunging region, we must choose a sensible boundary condition for the magnetic 
field on the marginally stable circular orbit. 
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The magnetic field on the marginally stable circular orbit is determined by the MHD 
processes in the disk region. In the disk region (r > r ms ), particles move on nearly circular 
orbits with \vf\ <C \vi\. As a disk particle moves a finite distance in the radial direction, 
it has undergone an infinite number (^> 1) of turns around the central black hole. As a 
result, the magnetic field lines frozen to the fluid are wound up around the center an infinite 
number of times. Therefore, in a stationary state, we expect that the magnetic field in the 
disk region is predominantly toroidal. In the Appendix we solve Maxwell's equations in the 
disk region and show that in the stationary state the magnetic field in the disk is likely 
to be parallel to the velocity field: B^/B r = u^/u r . Thus, in the stationary state, on the 
marginally stable circular orbit the magnetic field does not take any orientation but the one 
that satisfies the following condition 

1 . (82) 



On the marginally stable circular orbit, Eq. (^) implies that \Bf\ <C \B^\ since \vf\ <C 
From Eqs. (|48aj) and (|48b|) , we have 



-—7 = 1 + = r-jr . 83 

B r u<t> (-C u t + ^lu^)u* 

Then, Eq. (|52]) is satisfied if and only if \& = on the marginally stable circular orbit. In 
the stationary state $ is a constant, so we have 

* = (84) 



throughout the plunging region. Thus, in the stationary state, Eq. fl32|) holds throughout 
the plunging region 

Inserting Eqs. (^) and ( |3"3"D into Eq. (B2|), we obtain 

Bf Erv f 

B } TAV\n-uv})' { ' 

We have calculated 771, r/ 2 , and rj 3 in the plunging region, assuming that they all are 1 
at r = r ms and the boundary condition fl84D is satisfied. At r = r ms the toroidal magnetic 
field is assumed to be B^ = 0.05 in units of (47r/) m ) 1//2 , the corresponding radial magnetic 
field Bf is then given by Eq. (|85"D. The parameter 5 is taken to be 10~ 2 , 10~ 3 , 10~ 4 and 



10 , alternatively. The results are shown in Figs. 10 and O. From these figures we see 



that, in the plunging region the evolution of 771, r] 2 and r] 3 sensitively depends on the value of 
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5. For very small 5, T]i, 772 and 773 quickly decrease as the fluid gets into the plunging region. 
This is caused by the fact that a very small 5 (i.e., a very small \vfo\ oc 5 l l 2 ) corresponds 
to a very small \B? \ according to the boundary condition (p2[) , while the magnetic field in 
the plunging region is predominantly determined by \B?o\ instead of |-Bi | (see Figs. |3| and 
^). For a moderate 5 (e.g., 5 = 10 -2 ), rji, 772 and 773 may increase in the plunging region. 
However, even in this case, the conditions in Eq. ( [75| ) are always satisfied throughout the 
plunging region. 



B. Self-consistent solutions to the dynamical equations 

If we insert the solutions of Maxwell's equations that we obtained in Sec. [TV] into Eqs. ( |75|) 
and flTSI), then apply Eq. (|T7j), we obtain 

± (F m L + Anr 2 T EM /) = , (86) 

^ (F m E - 4nr 2 T EM /) = 0, (87) 

where 

w = ^ + ^ ' Tem '* r = §- T ^ r ■ m 

When \1/ = [i.e., Eq. ( jS2|) is satisfied in the plunging region], we have T EMt r = 0. Then 
Eq. (|87D implies that E keeps constant in the plunging region. However, since T EM J' 7^ 0, by 
Eq. (|86"D L varies in the plunging region. Therefore, for the solutions satisfying the boundary 
condition (|^) [or, equivalently, Eq. (|84])], the specific energy of particles is conserved but 
the specific angular momentum changes. 

Setting L = u^, E = —u t , and \1/ = 0, Eqs. (|86| ) and (^) can be integrated to obtain 

«* + ^ = -A, (89) 



u t = /e , (90) 



where 



co = - 7 — , h = =pr> &=Z]F' (91) 

V 1 m A m - 1 m 

where F^ = 47rr 2 J£ and F E = 47rr 2 J E are constants measuring the angular momentum flux 
and the energy flux across radius r, respectively. Using = g^u^ + g^u t , we can solve 

23 



Eqs. flU) and © for U<t> 



2Mac\ r r 

l»3^r JE JL 

1 + (1 _ 2M) 



(92) 



Then, from g u a Ub = — 1 and u r = g rr u r , we obtain the equation for u T 



\ru 



-AM 



1 - 



2A/ 



1 - 



2M 



2Ma 



1 + ^(1-^) . 



(93) 



Though the factor (l - ^M) appears in the denominators on the right-hand side of 
Eq. (|9~3"D, r = 2M is not a singularity of the equation, since the factor (l — —J disappears 
from the denominators if we expand the second term on the right-hand side of Eq. (|93l), then 



combine with the first term. However, the factor 



1 + (l — — ) 



of Eq. 



at 



r 2 u r 



2M 



. 



represents a singularity 



(94) 



The differentiation of Eq. ( p3|) with respect to u r gives rise to another singularity of the 
equation, which is at 



u r u — 



r 2 u r 



1 



2M 



This singularity appears as one differentiates Eq. ([ 
equation for u r . 

If we define the relativistic Alfven velocity by 



fi)=0- (95) 
) with respect to r to obtain a differential 



B 2 



(96) 



then we have 



1 + 



r 2 u r 



2M 



1 - r 2 



CAt-Ca 

u r u r 



and 



r 2 u r 



r 



„2 ' 



(97) 



(98) 



where c\ = CA a c A a < 1. Therefore, the singularities given by Eqs. (0) and (|95|) are critical 
points related to the Alfven speed: the Alfven point [Eq. ([94D 1 and the fast critical point 
[Eq. (§]i|. 



24 



Since we have assumed that the gas pressure is zero, the slow critical point is at u r = 
which is not relevant to us here. The Alfven point is not an X-type singularity and it 
does not impose any additional conditions on the solution for u r except setting the integral 
constant P5| |51fl . Therefore, what is really relevant here is the fast critical point given by 



Eq. (P3|). For weak fields, we expect that the fast critical point is located at a radius close to 
the marginally stable radius, where the accretion flow transits from subsonic motion in the 
disk region to supersonic motion in the plunging region || [7], ^J. Any physical solutions must 
smoothly pass the critical points, which sets strict constraints on the integral constants. 

It is far beyond the scope of the current paper to fully explore the properties of critical 
points in detail. However, as a first order approximation, we can assume that rf pa r ms , 
where rf is the radius at the fast critical point. Then, from Eq. (|95|) we have the radial 
velocity of the fluid at r = r ms 



Wq pa 



4\ns ( 2M 2 

' ms V ' ms 



1/3 



(99) 



If we set r = r ms in Eq. (|9"5|), then substitute Eq. ( |9"9"D into Eq. (|9"B|), we can solve for the 
constant /l, as a function of Co and /e- Therefore, for solutions that smoothly pass the 
critical points, among the three integral constants c , /e, and /l, only two are independent. 

Without explicitly solving the radial flow equation, i.e. Eq. (|93l), we can also obtain some 
interesting results on the horizon of the black hole. Since A — ► as r — * th, from Eq. (p3|) 
we have 

r 2 n 2 9M 

v^ = ^j- + — (/e-CWl) , (100) 

where = u r (r = r^). Then, from Eq. (|92|), we can obtain the specific angular momentum 
for particles at r = th 

Lh = -/ L + 4 c o • (101) 

Eq. (^) says that, with the boundary condition given by Eq. fl5I|), the specific energy 
of fluid particles does not change in the plunging region. However, the specific angular 
momentum of fluid particles does change. The specific angular momentum for particles on 
the marginally stable circular orbit, L , can be calculated from Eq. fl92|) by setting r = r ms . 
Then, by comparing Lh to L , we can see how much has changed in the specific angular 
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momentum as particles move from the marginally stable circular orbit to the black hole 
horizon. We have calculated the ratio 

AL 1 / , r N , N 

— = — (Lh - A)) , (102) 



and presented the results in Fig. 12 



To make the solutions smoothly joined to those in the disk region, in our calculations we 
have adopted that the specific energy is equal to that of a particle moving on the marginally 
stable circular orbit. Then, we have /e = — E ms . From Eqs. (|98| ) and (p9|), on the marginally 
stable circular orbit there is a one-to-one correspondence between Cq and 

= = T^hr • (103) 

47T/? m 1 - 4 

Therefore, we use T]i at r = r ms to specify the boundary value of the magnetic field, then cq 
is determined (up to a sign — which is not relevant here). Then, with the above approach, 
/l is determined as a function of rji and /e- We have chosen r]i to be 10~ 2 , 10~ 3 , and 10~ 4 , 
alternatively. 

From Fig. [12] we see that, though L is not constant in the plunging region, its variation 
is extremely small, if on the marginally stable circular orbit the dynamical effects of the 
magnetic field on the motion of the fluid particles is negligible (then r\\ <C 1 at r = r ms ). 
As a/M increases from —1 to 1, the variation in the specific angular momentum decreases, 
caused by the fact that the coordinate distance between the marginally stable circular orbit 
and the event horizon decreases with increasing a/M. For a black hole with a/M > 0, we 
have \AL/L\ < 3% if rji < 1CT 2 on r = r ms . For the extreme case of a/M = 1, we have 
\AL/L\ = 0, i.e. the specific angular momentum does not change at all. For the extreme 
case of a/M = —1, which has the largest coordinate distance from the marginally stable 
circular orbit to the event horizon, the variation in L is largest. However, even in this case, 
the variation in the specific angular momentum is also not big: \AL/L\ < 17% if rji < 10~ 2 
on r = r ms . 

For the same models we have also calculated the ratios rji and c^ r c A r /u r u r at r = th, 



and presented the results in Fig. [13]. Because of Eq. ( |103|) , the value of 771 is given by the 



left-hand side of Eq. (]98l) . The ratio CA r c A r / u r u r can be calculated by 

(104) 



Ur^ C 2 (/|-l+ 2 f) 
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r 2 u r 



From Fig. [13] we see that both B 2 /4:7ip m (solid curves) and c^ r c A r ju T u r (dashed curves) are 



small at r = rn- For a/M > 0, both are < 0.03 if r\\ < 10 2 at r 



The ratios go down 



as a/M increases. Even for the extreme case of a/M = —1, the ratios are also not big at 



rn if r7i < 10 2 at r 



B 2 /47rp m < 0.18 and c Ar c A r /u r u r < 0.07. 



Eq. ( |100| ) implies that is always finite. Then, u r u r - 
since A(r = th) = 0. Then, Eq. ( |104p implies that CA r c A ' 
right-hand side is finite. However, from Eq. (p7|) we have 



(r 2 /A) (u r ) 2 — > oo at r = Th 
-> oo also at r = rn since its 



CArC A 



U r U' 



< 1 



r=r H 



since 1 — 2M/r H = — a 2 /r^ and «g < . 

Since u r R is finite, Eqs. (|9S|) and ( |103| ) imply that 7/1 is also finite at r 
u r u r — > oo at r = th, we must have 

- 4) 



(105) 



?"h. Then, since 







(106) 



r=r H 



Eqs. ( 105 ) and (|106| ) imply that fluid particles always supersonically fall into the black hole. 



The numerical results in Fig. ^ confirm this conclusion. 

From Figs. IU - |T^, we can demonstrate that when the boundary condition fl52"|) [or, 
equivalently, (|84l) l is satisfied, the dynamical effects of magnetic fields in the plunging region 
are unimportant if they are so on the marginally stable circular orbit. 



VII. SUMMARY AND DISCUSSION 



With a simple model we have studied the evolution of magnetic fields in the plunging 
region around a Kerr black hole. The model contains the following assumptions: (1) The 
background spacetime is described by the Kerr metric; (2) The plasma is perfectly conduct- 
ing so that the magnetic field is frozen to the plasma fluid; (3) The kinematic approximation 
|T0| applies, i.e. the dynamical effects of the magnetic field on the fluid motion are negligible 
so that the plasma fluid flows along timelike geodesies toward the central black hole; (4) In 
a small neighborhood of the equatorial plane (i.e., \tt/2 — 9\ 1) the magnetic field and the 
velocity field have only radial and azimuthal components. The assumption (4) implies that 
u e = B e = and du /d8 = dB /d9 = on the equatorial plane. With above assumptions, 
we have exactly solved Maxwell's equations for axisymmetric solutions (i.e. d/dcf) = 0). 
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The solutions are given by Eqs. ( |48a|) and ( |48b|) [or, equivalently, Eqs. ( |50a|) and (|50b| )1, 
where Cq is a constant measuring the magnetic flux through a circle in the equatorial plane, 
\1/ is a function determining the time-evolution of the magnetic field. Both Cq and \1/ are 
determined by the boundary conditions. 

The dependence of the solutions on the coordinate time t is given by the function \1/ = 
ty(t — t), where r = r(r) is defined by Eq. (fJ2j). The function \1/ is uniquely determined by 
the initial and boundary conditions of the magnetic field. The general form of \l/ determines 
the retarded nature of the solutions: at any time the state of the magnetic field at a radius 
in the plunging region is determined by the state of the magnetic field at a larger radius and 
an earlier time (Figs. |l]and §). This suggests that the stationary state of the magnetic field 
in the plunging region is uniquely determined by the boundary conditions at the marginally 
stable circular orbit. 

Examples for the evolution of magnetic fields in the plunging region are shown in Figs. |3| 
- H|, for both the stationary case = constant, Figs. |3] - |5|) and the non- stationary case 
(\1/ 7^ constant, Figs. || - |]). From these figures we see that, the boundary value of the radial 
component of the magnetic field at the marginally stable circular orbit is more important 
in determining the strength of the magnetic field in the plunging region than the toroidal 
component. The initially toroidal component is attenuated in the plunging region by the 
radial expansion of the fluid (Fig. |3|). The initially radial component is amplified by the 
azimuthal and vertical compression (the convergence of the fluid), and a toroidal component 
is generated from the radial component by the shear motion of the fluid though the toroidal 
component is initially zero at the marginally stable circular orbit (Fig. §). This leads to 
the amplification of the magnetic field in the plunging region (Fig. |^). The evolution of the 
magnetic field depends on the spinning state of the black hole, but is more sensitive to the 
initial value of the radial velocity on the marginally stable circular orbit [or, equivalently 
the parameter 5 defined by Eq. fl57|) ]. 

The time-evolution of magnetic fields shown in Figs. |6] - |8] confirms our claim that the 
stationary state of the magnetic field in the plunging region is uniquely determined by 
the boundary conditions on the marginally stable circular orbit. If in the plunging region 
a magnetic field is initially deviated from the stationary solutions, it will evolve to and 
finally get saturated at the state given by the stationary solutions in a dynamical time scale 
determined by the free-fall motion of fluid particles in the plunging region. If the magnetic 



28 



field on the marginally stable circular orbit is in a stationary state, the magnetic field in 
the plunging region will automatically settle into a stationary state. Thus, the evolution of 
magnetic fields in the plunging region is very different from that in the disk region, in the 
latter case the state of the magnetic field is determined by the local MHD processes. The 
difference is caused by the following fact: in the disk region the fluid has a very small radial 
velocity so that local MHD processes have shorter time scales than the radial motion; in 
the plunging region the fluid has a large radial velocity so that local MHD processes usually 
have longer time scales than the radial motion. 

In deriving the solutions we have assumed that the dynamical effects of the magnetic field 
in the plunging region are unimportant [the kinematic approximation adopted in assumption 
(3)]. To justify this assumption, we have studied the dynamical effects of magnetic fields 
in the stationary state in two ways. First, we estimate the dynamical effects of magnetic 
fields by considering the back-reaction: substituting the solutions of Maxwell's equations 
we obtained, where we assumed that fluid particles move on geodesies, into the dynamical 
equations to check if the motion of the fluid is significantly affected by magnetic fields. In 
this way, we estimate the dynamical effects of magnetic fields on the motion of the fluid by 
calculating the parameters rji, r] 2 and rj 3 defined by Eq. ([78]). The mass density of the fluid, 
which appears in the denominators of r]i, 772 and 773, drops quickly in the plunging region 
if \vf\ -C I i^l on the marginally stable circular orbit (Fig. ^). However, the evolution of 
the magnetic field, which appears in the numerators of rji, 772 and 773, sensitively depends 
on the orientation of the magnetic field on the marginally stable circular orbit, the latter 
is essentially determined by the MHD processes in the disk region. If we require that the 
solutions in the plunging region are smoothly joined to the solutions in the thin Keplerian 
disk region, then it is reasonable to assume that on the marginally stable circular orbit, as 
well as in the disk region, the magnetic field is parallel to the velocity field: B r '/B^ = u r /v? 
(see the Appendix). This implies that for the solutions in the plunging region we should 
have \& = [Eq. fl8"4|)1, i.e. the orientation of the magnetic field follows the orientation of the 
velocity field of the fluid [Eq. (|8~2"D1. With such a boundary condition, we have calculated 771, 
772 and 7/3 , assuming that particles move on geodesies. We see that, if they are -C 1 on the 
marginally stable circular orbit, then they remain so in the plunging region (Figs. [10] and 
|TT] ). Indeed, 771, 7/2 and 773 decrease in the plunging region for sufficiently small on the 
marginally stable circular orbit. 
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Second, we self-consistently solve the coupled Maxwell and dynamical equations on the 
horizon of the black hole, check how much has changed in the specific angular momentum 
of fluid particles since they leave the marginally stable circular orbit. The specific energy of 
fluid particles is not changed by a magnetic field that satisfies $ = [see Eq. fl9"U|)]. However, 
the specific angular momentum does change [see Eq. Q89D]. We have calculated the specific 
angular momentum of fluid particles as they reach the horizon of the black hole, assuming 
that the fluid particles pass the fast critical point near the marginally stable circular orbit. 
The deviation from the specific angular momentum as the particles just leave the marginally 



stable circular orbit is shown in Fig. |T2j. We see that, though the specific angular momentum 
of fluid particles is changed by the magnetic field, the effects are always small assuming that 
on the marginally stable circular orbit the dynamical effects of the magnetic field are not 
important. For example, for a/M > we have \AL/L\ < 3% if rji < 10~ 2 on the marginally 
stable circular orbit. We have also calculated the ratios r/i = B 2 /Anp^ and c^ r c£ ju r u r 
on the horizon (Fig. [13]), where c A r is the radial component of the Alfven velocity. These 
two ratios are relevant to the critical points in the flow [see Eqs. fl97|), (|98D, and (|103| )1 and 
measure at what level the magnetic field affects the motion of the fluid. Fig. ^ shows that 
their values on the horizon are quite small. 

Both approaches (back-reaction and self-consistent solutions) confirm that the dynamical 
effects of the magnetic field are unimportant in the plunging region if they are so on the 
marginally stable circular orbit. 



Our results differ from that of Gammie ||12j| , in which he claimed that in the plunging 
region the dynamical effects of magnetic fields can be important. The difference is caused by 
the fact that Gammie used a boundary condition that is different from ours for the magnetic 
field. Gammie assumed that ^ = — Q ms Co, where f2 ms is the angular velocity of the disk 
at the marginally stable circular orbit, while we assume that \P = 0. Gammie's boundary 
condition implies that v$ = at r = r ms (as clearly stated in his paper), which makes r = r ms 
a singular point where p m = oo (to keep the mass flux F m nonzero). Hence, Gammie's solu- 
tions are not well-behaved at the marginally stable circular orbit. Our boundary condition 
allows the solutions in the plunging region to be smoothly joined to the solutions in the 
disk region (see the Appendix), since in our solutions all physical quantities are finite at 
the marginally stable orbit. Certainly, to precisely take care of the transition from the disk 
region to the plunging region, gas pressure and non-electromagnetic stress must be properly 
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taken into account near the inner edge of the disk. 

As mentioned in the Introduction, recently the "no-torque inner boundary condition" for 
thin accretion disks has been challenged by some authors (including Gammie) based on their 
studies on the evolution of magnetic fields in the plunging region [fET] , |T2|, [T^, [T7|]. The results 
in this paper suggest that in the plunging region the dynamical effects of the magnetic field 
are not important, if the solutions in the plunging region are smoothly joined to the the 
solutions in the thin Keplerian disk region. Thus, the argument against the "no-torque inner 
boundary condition" is not founded. 

Finally, we note that in the paper we have neglected dissipative processes like magnetic 
reconnection and ohmic dissipation, and the evaporation effect arising from the magnetic 
buoyancy and MHD instabilities. These processes operate in disks to limit the amplifica- 



tion of magnetic fields [10, ^3, 36, and references therein]. Certainly it is conceivable that 
they can also operate in the plunging region to reduce the amplification effect of magnetic 
fields. If these processes are important, the results of this paper tend to overestimate the 
amplification of magnetic fields in the plunging region. Then, the dynamical effects of mag- 
netic fields should be weaker than that we have estimated without considering dissipative 
processes, which will strengthen our conclusion that the dynamical effects of magnetic fields 
are unimportant in the plunging region. 
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APPENDIX: MAGNETIC FIELDS IN THE DISK REGION 

In the disk region (r > r ms ), viscosity plays an important role in the dynamics of disk 
particles. The viscosity transports angular momentum outward and dissipates energy, which 
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leads to disk accretion. So, in the disk region, particles move on non-geodesic worldlines. 

In the vertical direction, the gravity of the black hole is balanced by the gradient of 
the total pressure (= gas pressure + radiation pressure + magnetic field pressure) in the 
disk. Therefore, in a small neighborhood of the disk central plane, disk particles are more 
likely moving on planes parallel to the equatorial plane, rather than moving radially as in 
the plunging region. To describe such a motion, cylindrical coordinates (t, r, 0, z) are more 
suitable, where t and <fi have the same meaning as those in the Boyer-Lindquist coordinates, 
to the first order r also has the same meaning as that in the Boyer-Lindquist coordinates, 
but z = r (f - 6) where |f - 6\ < 1 . 

In the cylindrical coordinates, in a small neighborhood of the disk central plane the Kerr 
metric can be written as @] 

ds 2 = -\l \dt 2 dtd(p + —dr 2 + — deft 2 + dz 2 , (A.l) 

\ r J r A r 

where A and A are defined by Eq. ( |PT|) with 9 = tt/2 . 

The Maxwell equations that we need to solve are again given by Eq. ([H]), but now we 
have x a = (t, r, <p, z) and a/— g = r. As in the case for the plunging region, we assume that 
the motion of fluid particles is stationary but the evolution of the magnetic field can depend 
on time, i.e., we let du a /dt = but keep dB a /dt in the Maxwell equations. Furthermore, 
we assume that in a small neighborhood of the disk central plane u z = B z = so that 
du z /dz = dB z /dz = on the equatorial plane. Finally, we adopt du a /d<p = dB a /d<p = 
because of the axisymmetry of the system. Then, on the equatorial plane Eq. (|3l|) is reduced 
to 

r r ( u *fl/» - uPb*)} + \r(u r BP - u l3 B r )] = . (A.2) 
at 1 or L J 

Eq. (|A.2| ) can be solved with the same approach as that used in Sec. |V|. The solutions 



arc 



B r = -(-C ,d«t + *d««) , (A.3a) 
r 

B* = J_[-C 0jdUtU * + (1 + u%)* d ] , (A.3b) 
ru r 



where Co,d is an integral constant, — ^d(^ — t) is the solution of Eq. (^0|). The corre- 
sponding B l is related to B r and B^ by Eq. (BH). 
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Eq. ( |A.3a|) can be written as 



* d =—trB r + C 0>d Ut). (A.4) 

Assume that as r — > oo the disk is Keplerian so that = L oc r 1 / 2 and u t = —E « 1, and 
ri?' r keeps finite, then we must have 

lim ^ d = . (A.5) 



Eq. (|A.5| ) has an important implication for stationary solutions. For stationary solutions we 



have \l/d = constant, then by Eq. (|A.5| ) we must have ^ d = 0. So, for stationary solutions 
we have 

B r = ^ , (A.6a) 
r 

ru r 

where we have used u t = —E. 

From Eqs. (|A~6a| ) and (|A~6bD we have B <l '/B r = u (t> /u r } i.e. in the stationary state the 
magnetic field lines are parallel to the fluid motion. 
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t 



FIG. 1: The retarded nature of the solution given by Eq. (44). A time variation in the magnetic 
field at radius ro will lead to the same variation in the magnetic field at radius n < vq at a later 
time, assuming the fluid moves toward smaller radii (i.e. toward the central black hole). The time 
delay is given by the time needed by a fluid particle moving from ro to r±. 



36 



A 



*2 (t-r) 



II 



?1 (T-t) 



FIG. 2: A sketch of the solution given by Eq. (flff), corresponding to the boundary conditions given 
by Eq. (^). The parameter r is defined by Eq. (|42]), which is the coordinate time needed by a 
particle moving from radius ro to radius r. So, r = corresponds to r = ro, r = oo corresponds 
to the black hole horizon r = th- The coordinate time is represented by t. In this diagram, the 
worldlines of fluid particles are represented by straight lines t = r + constant. If the boundary 
conditions are imposed on the axes as shown in the diagram [Eq. (|45|)], the solution in the region 
between the r-axis and the dashed line t = r (i.e., region I) is determined by the boundary condition 
on the r-axis, the solution in the region between the i-axis and the dashed line (i.e., region II) is 
determined by the boundary condition on the t-axis. 
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FIG. 3: Stationary evolution of an initially toroidal magnetic field in the plunging region around a 
Kerr black hole. Each panel corresponds to a different spinning state of the black hole, as indicated 
by the dimensionless parameter a/M. Each curve starts from the marginally stable circular orbit 
(r = r ms , right end), ends at the horizon of the black hole (r = rn, left end). At r = r ms , 
the magnetic field is purely toroidal and has a value B± , the fluid particles have specific energy 
identical to that of the marginally stable circular orbit, and specific angular momentum smaller 
than that of the marginally stable circular orbit by a tiny fraction 5 [Eq. (|57|)]. In each panel, 
each curve corresponds to different values of 5: 10~ 2 , 10 -3 , 10 -4 and 10 -5 (downward). The radial 
component of the magnetic field is always zero. 
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FIG. 4: Similar to Fig. | but the magnetic field is initially radial, i.e. Bf = Bf-o and B^ = at 
r = r ms- Solid lines represent the toroidal component of the magnetic field, dashed lines represent 
the radial component. Though initially the toroidal component of the magnetic field is zero, the 
shear motion of the fluid in the plunging region generates toroidal component from the radial 
component. The evolution of the radial component does not depend on the initial kinetic state 
of the fluid, so there is only one dashed line in each panel. The four solid lines representing the 
toroidal component of the magnetic field correspond to different values of 5: 10 -2 , 10 -3 , 10 -4 and 
1CT 5 (upward). 
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FIG. 5: Stationary evolution of B 2 = B a B a (solid curves). The boundary conditions are the same 
as that in Fig. ^. The thin dashed line in each panel shows the absolute value of the radial velocity 
(i.e. \vf\) of the fluid, corresponding to 5 = 10 -5 . If we choose a different value of 5, the curve for 
the radial velocity will change slightly: the right end will change according to \v?\ oc (5 1 / 2 , the left 
end always approaches 1 (i.e. the speed of light). (For brevity, lg [B 2 /Bf ] and \vf \ use the same 
scale as labeled on the left side of the box in each panel.) 
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FIG. 6: Time evolution of magnetic fields around a Schwarzschild black hole (i.e. a = 0). Each 
panel corresponds to a particular moment. Thick solid lines show the toroidal component of the 
magnetic field (i.e., B£), thick dashed lines show the radial component (i.e., B?). The thin (upper) 
lines show the corresponding stationary state solutions. The fluid particles have specific energy 
identical to that of the marginally stable circular orbit, and specific angular momentum smaller 
than that of the marginally stable circular orbit by a factor of S = 10 -3 . The first (left and up) 
panel shows the initial and boundary conditions of the magnetic field: on r = r ms (right end), 
Bi = 0, Bf. = BfQ for all time t > 0; in the plunging region (rn < r < r ms ), Bt = at t = [i.e., 
Eqs. ( 64a ) and (|64b| ) ; the corresponding Bf at t = is automatically determined by the solutions]. 
The figures show that, as time goes on, the magnetic field quickly approaches and saturates at the 
state given by the stationary solutions. The black dot in each panel shows the radial position of a 
fluid particle at each moment: initially the particle is at r = r rj 
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FIG. 7: Same as Fig. |6] but for a Kerr black hole with a/M = 0.99. 
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FIG. 9: Stationary evolution of the mass density in the plunging region around a Kerr black hole. 
Each panel corresponds to a different spinning state of the black hole, as indicated by the values 
of a/M. The right end of each curve corresponds to the marginally stable circular orbit (t — r sm ). 
The left end of each curve corresponds to the horizon of the black hole (r = rn). The kinetic 



boundary conditions are given by Eq. (57). In each panel the four curves correspond respectively 
to 5 = 10~ 2 , 1CP 3 , 10 -4 and 10 -5 (downward). The mass density is in units of p m o - the mass 
density at the marginally stable circular orbit. 
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FIG. 10: Stationary evolution of the ratio i]\ = B 2 /Airp m in the plunging region. Each panel 
corresponds to a different spinning state of the black hole, as indicated by the values of a/M. The 
kinetic boundary conditions are given by Eq. (p7|). The boundary conditions for the magnetic field 
are: on r = r ms we have Bi = 0.05 in units of (47r / o m o) 1 / 2 . The corresponding B? is calculated with 
Eq. (|85|) so that the magnetic field is always parallel to the velocity field. In each panel the four 
curves correspond respectively to 5 = 10 -2 , 10 -3 , 10 -4 and 10 -5 (downward). Each curve starts 
from the marginal stable circular orbit (right end) and ends at the horizon of the black hole (left 
end). 
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FIG. 11: Same as Fig. |To| but for the ratios 
\BtB r / [A-k p ra Eu r )\ (dashed curves). 
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FIG. 12: The change in the specific angular momentum of fluid particles as they approach the event 
horizon of the black hole [Eq. ( |102| )1 , as a function of the spin of the black hole. It is assumed that 
the fast critical point is at r ~ r ms . The magnetic field is parallel to the velocity field of the fluid, 
then the specific energy of fluid particles keeps constant in the plunging region, which we have 
assumed to be equal to the specific energy of a particle moving on the marginally stable circular 
orbit. The boundary value of the magnetic field at r = r ms is specified by the ratio r]\ = B 2 /A7rp m . 
The three curves corresponds to three different values of 771 at r = r ms : 10 -2 , 10~ 3 , and 10~ 4 
(downward) . 
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